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1 Introduction 

Several years of efforfilifU has culminated in_a consistent effective field theory 
description of the nucleon nucleon interactional The ultimate goal of this ende- 
vour is to construct a framework with which to systematically describe bound 
and unbound multi-nucleon systems as well as elastic and inelastic processes. 
This effort was initiated by Weinberg's pioneering work on the subject d where 
he proposed a power-counting scheme for local-operators involving two or more 
nucleons and the inclusion of pions. However, it was shown that Weinberg's 
power counting scheme is not consistent. Recently, a consistent power-counting 
scheme has been proposecO to describe NN scattering and applied to observ- 
ables of the deuteron. I will focus on these subjects in this talk and give an 
indication of what this program of study might achcive in the near future. 

2 Why Effective Field Theory? 

Effective field theory is a very powerful technique for dealing with systems that 
possess widely seperated length scales. In a system with just two length scales 
li and I2 (as an example), the ratio Q = can be used as a small expansion 
parameter. Usually systems possess more symmetries when Q = and for 
small Q a perturbative expansion exists that makes use of the symmetries of 
the unperturbed system. 

In the theory of strong interactions there is one intrinsic length scale, 
Aqcd- For processes that occur at distance scales that are small compared to 
Aqcd the appropriate degrees of freedom are the quarks and gluons with the 
QCD Lagrange density used to compute processes as a perturbative expansion 
in the strong interaction coupling constant, a s (/x). In addition, there is power 
series expansion in forming the matrix elements of the quark-gluon operators, 
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with an expansion parameter Aqcd/Q- It is often useful to impose the con- 
straints of chiral symmetry, arising from the smallness of the light quark masses 
compared to Aqcd, giving yet a third expansion parameter, m q j Aqcd- An 
explicit example is the semileptonic decay of b-flavored hadrons. The inclusive 
decay rate for the decay of a f?-meson to a charmed final state isE 
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The quantities A, Ai and A2 are nonperturbative matrix elements correspond- 
ing to the energy of the light degrees of freedom of the £?-meson, the 6-quark 
fermi-motion and the chromomagnetic interaction between the 6-quark and 
the light degrees of freedom, respectively. In contrast, for processes at inter- 
mediate length scales there are no small expansion parameters. The strong 
interaction coupling a s ([i) is approaching unity, as is Aqcd/Q- At long dis- 
tances the appropriate degrees of freedom are the hadrons themselves, and not 
the quark and gluon fields. C/(l) em gauge symmetry of electromagnetism, the 
SU(2)l (g) SU(2)k chiral symmetry and the small expansions parameters p/A 
and m q /A (where p denote the external momentum of the hadrons, m q are the 
light quark masses and A is the scale of strong interactions) can be encorpo- 
rated a Lagrange density describing the low energy dynamics of hadrons and 
photons. 

Writing the Lagrange density as 

C(1^)=Y. C{1) ^ 0{1) ^ • ( 2 ) 

i 

where the operators (/x) are renormalized at the scale \i and are constructed 
from the hadronic fields. Coefficients Cw(/x) renormalized at jj, are determined 
by the short- distance behavior of the strong interactions. Clearly eq.(||) rep- 
resents an explicit seperation of scales. By construction, observables do not 
depend upon the scale /j, at which one chooses to renormalize. Initially, this 
appears to be a disasterous situation, since there are an infinite number of 
operators that one can construct and hence there are an infinite number of 
constants unconstrained by the symmetries of the theory. However, we have 
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learned much about dealing with non-renormalizable field theories over the last 
many years and in general such theories are predictive. To construct relations 
between certain observables that are valid to a certain precision, only a finite 
number of constants in the effective Lagrange density need to be determined. 
The key to making such theories predictive is to establish a consistent power 
counting scheme, one in which terms that "look small" at the level of the la- 
grange density do in fact make small contributions to observables, even with 
the inclusion of loop effects. 



3 The NN Interaction and Weinberg's Power Counting 

The most general Lagrange density consistent with chiral symmetry describing 
the interaction of two nucleons is 



f 2 f 2 
C = N^(iD + D 2 /2M)N + ^-Tr d^E^E + ^-wTrm,(E + E f ) 

v 4 ' 



- \c S {N^N) 2 -^C T {N^Nf+g A N^A-uN + ... , (3) 

where denotes a chiral-covariant derivative and E is the exponential of the 
isotriplet of pions 

E = exp (jM 
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with / = 132 MeV the pion decay constant and A is the axialvector meson field. 
The ellipses denote terms with more spatial derivatives and more insertions of 
the light quark mass matrix, m q . The Georgi-Manoha£3 naive dimensional 
analysis arising from a consideration of loop contributions to observables sug- 
gests that Cs,t ~ l// 2 , and Weinberg's power-counting will follow directly. 

A necessary ingredient for an EFT is a power counting scheme that dictates 
which graphs to compute in order to determine an observable to a desired 
order in the expansion. We denote the expansion parameters of the theory by 

1/2 

Q ~ |p|,m q . The main complication in the theory of nucleons and pions 
is the fact that a nucleon propagator S(q) = i/(qo ~ q 2 /2Af) scales like 1/Q 
if qo scales like m, or an external 3-momentum, while S(q) ~ M/Q if qo 
scales like an external kinetic energy. Similarly, in loops J dqo can scale like 
Q or Q 2 /M, depending on which pole is picked up. To distinguish between 
these two scaling properties it is convenient to define generalized "n-nucleon 
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Figure 1: The first two terms in the EFT expansion of the Feynman amplitude for nucleon- 
nucleon scattering in Weinberg's power-counting scheme. The leading amplitude Aq consists 

(2) 

of the sum of ladder diagrams with the leading 2-nucleon potential Vq at every rung. 

potentials" comprised of those parts of connected Feynman diagrams with 
In external nucleon lines that have no powers of M in their scaling (except 
from relativistic corrections). y(") includes diagrams which are n-nucleon 
irreducible and parts of diagrams which are 1-nucleon irreducible. To compute 
the latter contribution to K'"' one identifies all combinations of two or more 
internal nucleon lines that can be simultaneously on-shell, and excludes their 
pole contributions when performing the J dqo loop integrations. 

Two nucleon scattering is simple since the graphs are all ladder diagrams 
with insertions of yw's acting as ladder rungs. Each loop of the ladder in- 
troduces a loop integration (dgod 3 q ~ Q 5 /M) and two nucleon propagators 
(M 2 /Q 4 ) to give a factor of (QM) per loop. If one treats M ~ Q°, it follows 
that perturbation theory is adequate for describing the 2-nucleon system at 
low Q. In order to accommodate large scattering lengths and bound states 
near threshold, as in the ^So and 3 5i — 3 Di channels one must conclude that 
M ~ l/Q in this power counting scheme. At leading order one must sum up all 

(2) 

ladder diagrams with insertions of the leading two-body potential V Q , while 

(2) 

at subleading order one includes one insertion of the subleading potential, Vy ; 

(2) 

and all powers of Vq , and so forth. 

At leading order in Weinbergs power-counting there are contributions to 

(2) 

Vq from both the local four-nucleon operators, Cs,t and from the exchange 
of a single potential pion, giving a momentum space potential of 

V {2) (t> t>')-C ( 9 * a \ (q-^q-^Xn-ra) f . 

(P ' P) " C WJ (* 2 +™l) ' (5) 

where C denotes the combination of Cs,t appropriate for a given spin-isospin 
channel. The leading order amplitude results from summing the graphs shown 
in Fig. ([[]) resulting from this potential, i.e. solving the Schrodinger equation. 
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(a) 



(b) 



Figure 2: Graphs with logarithmic divergences. The divergence in graph (a) is proportional 
to M 2 m 2 , while graph (b) has a divergence proportional to Af 2 p 4 . The solid lines are 
nucleons and the dashed lines are pions. 



Figure 3: A contribution to the pion ladder sum, arising at leading order in Weinberg's 
power-counting. The solid lines are nucleons and the dashed lines are pions. 



In the So channel at two-loops in the ladder sum there is a logarithmic diver- 
gence in the graph shown in Fig. (^a) that must be regulated. In dimensional 
regularization the divergent part of this graph is 

l g 2 A mlM 2 - 
"7T28^ C ' (6) 

which requires a counterterm with a single insertion of the light quark mass 
matrix. However, the coefficients of these operators must scale like M 2 , and 
since m^M 2 ~ Q° these formally higher order operators in Weinberg's power- 
counting are required at leading order to absorb divergences in the time-ordered 

(2) 

products of the leading order potential, Vq . Ignoring the multi-pion vertices 
arising from these operators, they can be re-absorbed into the leading operators 
with coefficients Cs,t- 

The situation is different in the 3 S*i — 3 Di channel and in higher partial 
waves. A contribution to the leading order ladder sum is shown in Fig. (^|), 
arising from seven potential pion exchanges, i.e a six-loop graph. It is straight- 
forward to deduce that this graph has a logarithmic divergence at order (QM ) 6 , 
and therefore, counterterms involving V 6 are required at leading order in the 
expansion. Clearly, the same discussion can be made for an arbitrary number 
of potential pion exchanges, and therefore counterterms involving an arbitrary 
even number of V's are required. This is a clear demonstration of the failure 
of Weinbergs power-counting. Further, this conclusion is true for all regular- 
ization schemes and not just for dimensional regularization. 



5 



4 A New Power Counting 

Lets us begin by examining the general form of the amplitude for nucleon 
scattering in a S-wave 

A= ^L I . (7) 

M p cot S — ip 

From quantum mechanics it is well known that p cot 5 has a momentum ex- 
pansion for p <C A (the effective range expansion) , 

1 1 °° / r?\ n+1 
pcoa= __ + A 2 ]Tr„ , (8) 

n=0 ^ ' 



where a is the scattering length, and tq is the effective range. For scattering in 
the ^So and 3 Si channels the scattering lengths are found to be large, = 
-23.714 ± 0.013 fm and a^ 51 ' = +5.425 ± 0.0014 fm respectively. Expanding 
the expression for the amplitude in eq.(^) in powers of p/A while retaining ap 
to all orders gives 



4tt 1 

A = - 



M (1/a + ip) 



ro/2 v + /; o/2) ; 9 / + /;; /2A \/ 



(l/a + ipy (1/a + ip) 2 ^ (1/a + ip) 

(9) 

For p > l/\a\ the terms in this expansion scale as {p . . .}, and the 

expansion in the effective theory takes the form 

oo 

A= ]T A n , A n ~ p n . (10) 
n=-l 

In the theory without pions we can explicitly compute the s-wave ampli- 
tude in each spin channel to all orders in the momentum expansion, 

a = £ ° 2np2n fm 

1 + M(/i + zp)/4^£ C 2nP 2n ' 1 1 

where Cm is the coefficient of the p 2n term in the lagrange density, /i is the 
mnormalization scale and we have used Power Divergence Subtraction (PDS 
p\ to define the theory. A typical loop graph that appears in the amplitude 
has the form 



= -M(MEn-ME-ur- 3)/2 r (^) ( ip-i)/ 2 • ( 12 ) 
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Figure 4: Leading and subleading contributions arising from local operators. 

where D is the number of space-time dimensions. In the PDS scheme the pole 
at D = 3 is removed by adding a local counterterm to the lagrange density, so 
that the sum of the loop graph and counterterm in D — 4 dimensions is 

In DS = In + SI n = -(ME)* f^) (p + ip). (13) 

The amplitude A is independent of the subtraction point fi and this deter- 
mines the n dependence of the coefficients, C2n- In the PDS scheme one finds 
that for /i 3> l/|a|, the couplings C l 2n(A t ) scale as 

4-7T 

CM - j^j^ , (14) 

so that if we take fj, ~ p, C2 n (/i) ~ l/p n+1 . A factor of V 2 ™ at a vertex 
scales as p 2n , while each loop contributes a factor of p. Therefore, the leading 
order contribution to the scattering amplitude A-i scales as p _1 and consists 
of the sum of bubble diagrams with Co vertices. Contributions scaling as 
higher powers of p come from perturbative insertions of derivative interactions, 
dressed to all orders by Co- The first two terms in the expansion 

A _ -Co a -C2P 2 , , 

[1 + ^0* + *)] ' ° [l + ^ + ^p)] 2 ' ( } 

correspond to the Feynman diagrams in Fig. ([|). A comparison with eq.([|) 
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Figure 5: Contributions to the /3-functions for Co in the theory with pions 
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The dependence of C 2n (p) on is determined by requiring the amplitude be 
independent of the renormalization scale /x. The physical parameters a, r n 
enter as boundary conditions on the resulting renormalization group (RG) 
equations. In general one finds the coefficients to be 



-in 



M(-fi + I /a) 1/a 



ro/2 



0( M -) 



(17) 



which has the scaling property in eq.(14). The leading behavior depends on 
the two parameters a and ro encountered when solving for Cq([i) and C2(/i). 
This is due to the C 2n couplings being driven primarily by lower dimensional 
interactions. 

The inclusion of pions into the theory is straightforward. While the co- 
efficients of the local operators are renormalized, and scale as powers of the 



renormalization scale \i (we use Q = fi ~ p 



1/2- 



, the exchange of a single 



potential pion does not suffer from such renormalizations and therefore pion 
exchange is a sub-leading contribution, Q°. At the same order as the exchange 
of a potential pion is an insertion of a C 2 operator and a single insertion of the 
quark mass matrix m q . Ignoring isospin violation, these operators involving 
insertions of the light quark mass matrix with coefficients D 2 have the same 
structure as the Co operators. For \i ~ m„, Co(fi) cx C 2 {[i) cx 1/^i 2 and 
D 2 (fi) cx 1/fJ? for the 1 jSo and 3 S\ channels. A feature of the theory with 
pions is that this scaling behavior breaks down at low momentum, p ~ 1/H, 
and at sufficiently high momentum. Solving the RG equation arising from 
the graphs shown in Fig. (|5|) in the 1 So channel with the boundary condi- 

tion Cq (0) = 4nai/M, where a\ is the So scattering length, we find for 
A 4 > l/|«i| 



C, 
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Figure 6: Graphs contributing to the subleading amplitude Ao- The shaded ovals are defined 
in Fig. (Q). 



with 

Ann = ^ ~ 300 MeV , (19) 

and therefore the power counting changes when fi <~ Ann- The UV fixed point 
toward which Cq is driven largely cancels the ^-function component of the 
single potential pion exchange in the Sq channel. As a result, this power 
counting is valid up to p ~ Ann and the power counting in both channels 
is expected to fail at momenta on the order of Ann- We conclude that the 
expansion parameter for this theory is ~ m^/ Ann ~ 5, larger than one would 
like. 



5 NN scattering in the x 5o Channel 

Having established a consistent power-counting in the previous sections we now 
apply it to NN scattering in the ^So channel. The amplitude at order Q -1 
and Q° determined from the graphs shown in Fig. (0) and Fig. (0) are 
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Figure 7: The phase shift <5 for the 1 5o channel. The dot-dashed curve is the one parameter fit 
at order Q~ 1 , that reproduces the scattering length. The dashed curve corresponds to fitting 
8 between < p < 200 MeV, while the dotted curve corresponds to fitting the scattering 
length and effective range. The solid line shows the results of the Nijmegen partial wave 
analysis. 
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/la \ 

At order Q~ there is one unknown coefficient Cq d that must be determined 

/Xq \ /Xq \ 

from data while at order Q° there are three unknown coefficients Cq , 



and Z?2 that must be determined. The graph giving ^4q J v ' is divergent in 
four dimensions and therefore gives rise to the logarithmic dependence on the 
renormalization scale fx in eq. (po|). In order for the expansion to converge, the 
leading term A—\ must capture most of the scattering length. The phase shift 
S is perturbatively expanded in Q, S = 5^ + d^> + . . .. and fit to the results of 



(IV) 



10 




p (MeV) 




Figure 8: The phase shifts <5o, S2 and mixing parameter e\ for the 3 Si — 3 Di channel. The 
solid line denotes the results of the Nijmegen partial wave analysis. The dot-dashed curve 
is the fit at order Q~ 1 for So, while 62 = £1 = at this order. The dashed curves are 
the results of the order Q° fit of So to the partial wave analysis over the momentum range 
p < 200 MeV. 

the Nijmegan partial-wave analysisS over a momentum range p < 200 MeV. 
We find for /i = m n 

C*^ So) = -3.34 fm 2 , D ( 2 So) = -0.42 fm 4 , cf So) = 3.24 fm 4 , (21) 

giving the dashed curve plotted in Fig. (Q). It is clear from Fig. (Q) that the 
corrections to the leading order result become substantial above ~ 200 MeV 
and we expect the expansion to become unreliable at momenta larger than 
this value. We chose to renormalize at fj, = m v for our numerical analysis, but 
we could have chosen any value of with Ann S> /i ^> I /a. The logarithm 
appearing in the subleading amplitude suggests we choose /1 ~ . 

6 NN scattering in the 3 5i - 3 L>i Channel 

The analysis of scattering in the 3 Si— 3 D± channel is a straightforward extension 
of the analysis in the ^0 channel. The important difference is that the nucleons 
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in the initial and final states with total angular momentum J = 1 can be in an 
orbital angular momentum state of either L = or L = 2. The power counting 
for amplitudes that take the nucleons from a 3 Si-state to a 3 5i-state is identical 
to the analysis in the ^o-channel. Operators between two 3 Di states are not 
directly renormalized by the leading operators, which project out only 3 Si 
states. However, they are renormalized by operators that mix the 3 <Si and 3 D\ 
states, which in turn are renormalized by the leading interactions. Further, 
they involve a total of four spatial derivatives, two on the incoming nucleons, 
and two on the out-going nucleons. Therefore, such operators contribute at 
order Q 3 , and can be neglected in the present computation. Consequently, 
amplitudes for scattering from an 3 D\ state into an 3 D\ state are dominated by 
single potential pion exchange at order Q°. Operators connecting 3 D\ and 3 Si 
states arc renormalized by the leading operators, but only on the L — "side" 

of the operator. Therefore the coefficient of this operator, C\ 1 ' ~ 
contributing at order Q 1 and it can be neglected at order Q°. Thus, mixing 
between 3 D\ and 3 Si states is dominated by single potential pion exchange 
dressed by a bubble chain of Cq operators and a parameter free prediction 

fig \ /3g \ 

for this mixing exists at order Q°. Fitting the parameters Cq 1 , C% and 

£>2 to the phase shift So over the momentum range p < 200 McV yields, at 
/! = m n 

C^ Sl) = -5.51 fm 2 , L>f Sl) = 1.32 fm 4 , cf Sl) = 9.91 fm 4 . (22) 

The dashed curves in Fig. (@) show the phase shifts &i, 62 and mixing parameter 
£i compared to the Nijmegen partial wave analysisc 3 ] for this set of coefficients. 

( 3 s ) 

There are no free parameters at this order in either £1 or 62 once C 01 ; has 
been determined from Sq. 



7 The Deuteron 

Once the Lagrange density in the nucleoli sector has been established the 
standard-tools of field theory can be used to determine the properties of the 
deuterortZI. To compute the electromagnetic form factors of the deuteron one 
first computes the three point correlation function between a source that cre- 
ates a nucleon pair in a 3 jSi state, a source that destroys a nucleon pair in a 3 <Si 
state and a source that creates a photon. After LSZ reduction and wavefunc- 
tion renormalization one obtains the electromagnetic form factors. Leading, 
subleading and subsubleading order graphs contributing to the electric form 
factors of the deuteron are shown in Fig. (^). The resulting form factors A(q 2 ) 
and B(q 2 ) that appear in the differential cross section for electron deuteron 
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Figure 9: The diagrams contributing to the electric form factors of the deuteron. 



scattering are shown in Fig. (p0[). A(q 2 ) is dominated by the charge form fac- 
tor and B(q 2 ) depends only upon the magnetic form factor. One sees that 
the form factors computed at subleading order agree well with the data. The 
charge radius of the deuteron is found to be 

_ 1 , r (,A M ^-^ 2 + g 2 A Mml(?,m« + IO7) 

V<r >th ~^ + ^ + 288V2^(m. + 2 7 )3 + " 
= 1.89 fm+ ... 

V< r 2 > cxpt = 1.963 fm , (23) 



where 7 = \J AI^B is the binding momentum of the deuteron. The result 
computed to subleading order agrees with the observed value. The quadrupolc 
moment vanishes at leading order in the expansion but receives a contribution 
at subleading order from the exchange of one potential pion, giving 

g 2 A M(QMB + 9m w j + 4m 2 ) 



^ Q ' th 307r/2(m w + 2 7 )3 

= 0.40 fm 2 + ... 
IJ,Q,expt = 0.2859 fm 2 , (24) 

which is approximately 30% larger than the experimental value. Clearly a 
subsubleading calculation is needed to ensure that the theoretical value ia.con- 
verging to the experimental value. This work is currently in progres£3. In 
contrast, there is a local counterterm contributing to the magnetic moment at 
subleading order. We find 

MM = + f*n + L 2 ((J,) — (H ~ if + ... 

= 0.88 - 0.02 (fit) , (25) 
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Figure 10: The form factors A(q 2 ) and B{q 2 ) measured in elastic electron-deuteron scat- 
tering. The dashed curve is the leading order prediction, while the solid curve subleading 
prediction. There is one counterterm in B(q 2 ) that is fixed by the deuteron magnetic mo- 
ment. 



which determines the counterterm L 2 at the scale fi. Once this counterterm is 
determined from the deuteron magnetic moment, the entire form factor B(q 2 ) 
is determined to subleading order. 

We have also computed the polarizabilities of the deuteroiH, and find that 
the scalar cxeo and tensor oe2 electric polarizabilities, are 



«E0 




«£0 



OLE1 



+ 



ag\M% ml (3ml + Wm^j + 24 7 2 ) 
384tt/ 2 ^(m^ + 2 7 ) 4 + "' 

ag 2 A M 2 N 2ml + jj^47 + 16m^ 7 2 + 8 7 3 
80tt/ 2 7 2 (to w + 2 7 ) 4 

0.595 fm 3 + ... , a E 2 = -0.062 fm 3 + ... 



(26) 
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while the scalar and tensor magnetic polarizabiltities, (3 mo and [3m2, are 



a 



2M 



N 



16 7 2 



3 7 2 



Gtt 



0. 



M2 



Mjv 2 7 2 + 4tt 7 ^ 1 j 



P M0 = 0.067 fm d 



J/2 



0.195 fm' 5 



(27) 



where A_ ± a (E) is the scattering amplitude in the So channel evaluated at a 
center of mass energy E. and 

«;W are the isoscalar and isovector magnetic 
moments of the nucleon repectively. 



8 Present Limitations 

The power counting fails at the scale Ann and there has been little progress 
in understanding how to deal momentum higher than Ann in the effective 
field theory framework. Therefore we are presently unable to address some 
important areas, such as photo-pion production off the deuteron or elastic 
pion deuteron scattering. The typical momentum scale in such a process is 
~ s/Mp^rn^ greater than Ann- 

It is also worth commenting on the role of baryonic resonances in this 
program. The impact of the A resonance has been determined in Weinberg's 
power-counting scheme in two different prescriptionsrcl |^. It is found not to 
play an important role in NN scattering as the mass scale that sets the size of 
its contribution is ^/A/(A/a — M) ~ 500 MeV. This scale is higher than the 
scale at which the power counting becomes inappropriate, Ann, and so the 
baryonic resonances should not be included, until the theory above the scale 
Ann is constructed. 



9 Conclusions 

After several years of investigation we have constructed a consistent power 
counting for an effective field theory description of the nucleon nucleon interactio 
Pions are subleading compared to the local momentum independent four- 
nucleon operators and can be treated in perturbation theory. NN scattering 
in the 1 Sq and 3 Si — 3 -Di channels to sub-leading order has been considered 
and most impressive perhaps is the parameter-free prediction of the 3 Si — 3 Z?i 
mixing parameter e\ which agrees reasonably well with the Nijmegen phase 

c The A(1232)— and other baryon resonances have been consistently included in the single 
nucleon sectotfl. 
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shift analysis. The properties of the deuteron bound state follow straightfor- 
wardly from this construction. I presented the electromagnetic form factors 
and polarizabilities. One of the nice features of having a field theory construc- 
tion is that there are many well known theorems. One such theorenflla is 
that contributions from operators that vanish by the equations of motion are 
of the form of higher dimension operators that do not vanish by the equations 
of motion and therefore such operators can be neglected. 

The future looks extremely promising for a systematic analysis of nuclear 
physics using effective field theory. The short term program will be to continue 
to examine the two-body systems in detail, working to higher orders in the ef- 
fective field theory expansion. In the long-term one hopes to make progress 
in three-bodjE-a and higher-body systems. Both lines of investigation are di- 
rectly relevent to the BLAST program, and I look forward to close cooperation 
between experimental and theoretical endevours in this area. 



I would like to thank the organizers Ricardo Alarcon and Richard Milncr 
for putting together such a stimulating meeting. This work is supported in 
part by Department of Energy Grant DE-FG03-97ER41014. 
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